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ABSTRACT 

We report on the discovery of an isolated, compact HII region in the Virgo cluster. The object is 
located in the diffuse outer halo of NGC 4388, or could possibly be in intracluster space. Star formation 
can thus take place far outside the main star forming regions of galaxies. This object is powered by a 
small starburst with an estimated mass of ~ 400 M© and age of ~ 3Myr. From a total sample of 17 
HII region candidates, the present rate of isolated star formation estimated in our Virgo field is small, 
^ 10“® MQarcmin”^ yr“^. However, this mode of star formation might have been more important at 
higher redshifts and be responsible for a fraction of the observed intracluster stars and total cluster 
metal production. This object is relevant also for distance determinations with the planetary nebula 
luminosity function from emission line surveys, for high-velocity clouds and the in situ origin of B stars 
in the Galactic halo, and for local enrichment of the intracluster gas by Type H supernovae. 

Subject headings: stars: formation - HII regions - galaxies: ISM - galaxies: star clusters - galaxies: 
abundances - intergalactic medium 


1. INTRODUCTION 

Stars are usually observed to form in galaxies, that is, 
in disks, dwarfs, and starbursts. In radio galaxies, star 
formation may be triggered by energetic jet outflows (e.g., 

Bicknell et al. 2000). The HII region we describe here 
shows that isolated star formation takes place in the dif¬ 
fuse outskirts of galaxies, at the boundary of, if not already 
in, Virgo intracluster space. 

In nearby galaxy clusters, a diffuse intracluster star com¬ 
ponent has been inferred from surface brightness measure¬ 
ments (Bernstein et al. 1995) and detection of individual 
stars (e.g., Arnaboldi et al. 1996, Ferguson et al. 1998, 

Feldmeier 2002). Its origin may be explained readily by 
dynamical processes acting on low-surface brightness disks 
and dwarfs, which unbind the stars from these galaxies 
(Moore et al. 1999). However, gas may be efficiently re¬ 
moved from infalling galaxies by ram pressure stripping 
(Quilis et al. 2000, Gavazzi et al. 2001), or may be tidally 
dissolved, or could fall into the cluster as pristine clouds. 

Some of this gas may form stars, which would also con¬ 
tribute to the diffuse component. We discuss here the 
current rate of such star formation from isolated compact 
HII regions (ICHHs) and from jet induced star formation 
in a Virgo cluster field. 

2. AN HII REGION IN VIRGO: OBSERVED PROPERTIES 

The target object was found in an emission line survey 
for planetary nebulae in a Virgo intracluster field, centred 
at a(J2000)12 : 25 : 31.9, ,5(J2000)12 : 43 : 47.7, using 
Ho and [OHI] narrow band and V-l-R broad band filter 

^ Based on observations carried out at UT4 of the VLT, Paranal, Chile, which is operated by the European Southern Observatory. 

^No correction for atmospheric refraction is needed since FORS2 has an atmospheric dispersion corrector. 
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photometry (Arnaboldi et al. 2002, Okamura et al. 2002), 
with Suprime-Gam on the Subaru Telescope. Because of 
the large Ha to [OHI] flux ratio the (unresolved) object 
was classified as a candidate compact HII region. 

A spectrum was taken at UT4 of the VLT at Paranal, 
on the night of April 14, 2002, with FORS2 in MOS mode. 
The observations were carried out with GRISM-150I and 
the order separation filter GG435-I-81, giving a wavelength 
coverage of 4500 — 10200 A and a dispersion of 6.7 A pix“^. 
The slit width was 1.4 arcsec, and the angular scale along 
the slitlet was 0.126 arcsec pix“^. The total exposure time 
was 7 X 1800s. The nights were clear but not photometric; 
the mean seeing was better than 1.0 arcsec. Spectrophoto- 
metric standard stars were observed at the beginning and 
end of the night, but cirrus clouds at these times made the 
flux calibration uncertain. 

The data reduction was carried out using standard tasks 
in IRAF, using an arc lamp wavelength calibration and 
observations of a spectrophotometric standard star. The 
spectrum was corrected for atmospheric extinction, using 
a table for La Silla^. 

The wavelength and flux-calibrated spectrum is shown 
in Figure 1: it clearly shows a blue continuum and a 
number of emission lines at Virgo redshift: Ha, H/3, 
Hy, [OHI]A5007 and A4959, [SH]AA6717 -b 6731, and 
[SIH]A9069. We do not resolve Ha and [NHJA6548, and 
the [SH]AA6717,6731 lines, respectively, and we do not see 
the weaker [OI]A6300, [OHIJA4363, and [SHI]A6312 lines. 
The Hy line lies close to the blue end of the spectrum, and 
the [SHIJA9069 line in a region of strong sky emission. 
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Fig. 1.— Observed emission spectrum. 


Table 1 

Observed and reddening-corrected emission line fluxes relative to H/3. 


Line 

A 

A 

observed flux 
rel. to H/3 

de-reddened flux 
rel. to H/3 

H 7 

4340 

0.58 

0.71 

H/3 

4861 

1.00 

1.00 

[OIII] 

4959 

1.03 

0.99 

[OIII] 

5007 

3.05 

2.91 

[NI] 

5200 

0.08 


Ha 

6563 

4.30 

2.89 

[Nil] 

6583 

0.75 

0.50 

[SH] 

6717-b6731 

0.40 

0.26 

[SHI] 

9069 

0.30 

0.14 


[NI] identification uncertain. 

[NII]A6583 flux from two-Gaussian fit. 

Ha flux is corrected for corresponding [NIIJA6548 emission. 
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The observed emission lines and their uncorrected and 
reddening-corrected fluxes normalized to H/3 are listed in 
Table 1. The errors in the fluxes are approximately 10% 
of the H/3 flux, and larger at the blue edge of the spectral 
range. The Ha and [NIIJA6548 lines are unresolved, with 
the [Nil] A6583 line appearing in the red wing of Ha. From 
a two-Gaussian fit to the combined emission we estimate 
the line ratio [NH]A6583/(Ha-f [NH]A6548) ~ 0.17. This 
is consistent with the value expected for an HII region with 
the large observed [OHI]A5007/H/3 ratio; see below. The 
Ha flux given in the table is that measured as the total flux 
in the line reduced by one third that in [NH]A6583. The 
observed Balmer decrement Ha/H/3 then becomes 4.3. 

If the intrinsic Balmer decrement has the theoretical 
value of 2.89 for large optical depth (case B) and temper¬ 
ature T = 10^ K, then E(B-V)=0.40 and Av=1.23 mag. 
Most of this reddening is intrinsic to the source, as the 
galactic absorption to nearby galaxies in the Virgo cluster 
is about Ay — 0.12. Line fluxes were corrected for redden¬ 
ing using this value of E(B-V) and the reddening curve of 
Cardelli et al. (1989). Compared to H/3, the corrected Hy 
flux is about 1.5 times larger than expected; however, we 
cannot regard the Hy flux as reliable, because the line lies 
right at the edge of the spectrum where uncertainties in 
the grism response, sky subtraction and determination of 
the continuum level will be at their worst. 

The corrected line ratios [OHIJA5007/H/3 = 2.91, 
[NH]A6583/Ha = 0.17, and [SH] AA6717,6731/Ha = 0.09 
place the object clearly into the domain of HII regions 
in Figs. 1,2 of Veilleux & Osterbrock (1987; V087). 
This is also confirmed by the weakness (non-detection) of 
[OI]A6300 in the spectrum. 

Because the flux calibration using our standard stars is 
uncertain, we transformed it to the Jacoby (1989) mom 
calibration of the imaging photometry from Suprime-Cam 
(Arnaboldi et al. 2002). This gives a total mom = 25.7 
magnitude for this object. Because of the [OHI] fil¬ 
ter width used in the emission line survey, this flux in¬ 
cludes both the [OHI]AA4959,5007 lines, as well as the 
flux from the continuum. Using the spectrum to mea¬ 
sure the fraction of light in the continuum, we can correct 
to the magnitude in the A5007 line only: msooy = 26.15. 
From the standard star-calibrated spectrum, we would get 
1115007 = —2.51ogF5oo7 —13.74 = 25.1. The difference gives 
us a factor 2.63, by which to rescale F 5007 to match the 
photometric calibration. The rescaled total flux in the 
[OHI] A5007 A line is F 5007 = 11.0 x 10 ~^'^ergs“^ cm“^. 
With E(B-V)=0.4 the dereddened flux becomes ^5007 = 
4.0 X 10“^® ergs“^ cm“^. The fluxes in the other emis¬ 
sion lines can then be inferred from the dereddened line 
ratios in Table 1; in particular the Ha flux becomes 
3.9 X 10“^® ergs“^ cm“^. 

The measured V-band continuum flux from the spec¬ 
trum, when applying the same calibration and correcting 
for reddening, becomes Fy = 7.0 x 10“^® ergs“^ cm“^. 
The corresponding dereddened apparent magnitude is 
my = 24.2. Absolute magnitudes and fluxes will be com¬ 
puted for a nominal Virgo distance modulus of 31.16, cor¬ 
responding to 17Mpc distance (Tonry et al. 2001). How¬ 
ever, we will sometimes keep the distance dependence by 
writing D = dn Mpc. The continuum my = 24.2 then 
gives My = —7.0. 


3. PHYSICAL PARAMETERS 


Metallicity: The high values of [OHI]/Ha = 1.0 and 
[OHI]/H/3 = 2.9 indicate subsolar metallicity. The V087 
diagrams are not ideal for estimating metallicities, as dis¬ 
cussed by Dopita et al. (2000), but from Figs. 2,3 of that 
paper we estimate Z ~ 0.4 and an ionization parame¬ 
ter g ~ 4 X 10^ (logC/ = —2.9). Both very low and 
high metallicities are not consistent with their model re¬ 
sults and these line ratios. A more accurate determi¬ 
nation is in principle possible using the S '23 parameter, 
defined as S '23 =([SH]AA6717,6731-f [SHI]AA9069,9532) / 
H/3 = 0.75, where we have used the theoretical ratio 
[SHI]A9532/[SHI]A9069 = 2.48. From the empirical cal¬ 
ibration of Diaz & Perez-Montero (2000) we then find 
12-|-log(0/H) ~ 8.08±0.2 (about 0.15 to 0.25 solar with or 
without the depletion factor used by Dopita et al. 2000). 
In our case, S 23 has additional uncertainties due to the 
weak [SH] line and the bright NIR sky lines, which might 
affect the [SHI] flux. In the following, we therefore adopt 
Z ~ 0.4. 

Temperature, density, ionization parameter: Un¬ 
fortunately, the S/N of the spectrum is not large enough 
to detect the weak [OHI]A4363 and [NH]A5755 lines used 
to determine electron temperature, nor do we have lines 
to determine the electron density. In the following we use 
Te = 10^ K when needed. 

Estimates of the ionization parameter from the Sulfur 
lines, using equations (8) and (6) of Diaz et al. (2000) 
give similar values to that from the V087 diagrams: 
[SH]AA6717,6731/[SHI]AA9069,9532 ~ 0.53 gives logU = 
—2.5, and [SH]AA6'717,6731/H/3 ~ 0.26 and metallicity 0.4 
solar gives logU = —2.9. 

Luminosities: From above, the total V-band luminos¬ 
ity is 

Ly = dTrU^Uy = 2.4 X lO^^d^yergs-^ (1) 
The Ha flux is 

Lna = = 1.3 X 10®^di7 ergs"^ (2) 

The total number of H-ionizing photons (Osterbrock 1989) 
with recombination coefficients as and a|[^ for T = 10^ K 
is 

g(H°) = 2.m Lna/hiyua = 1-3 x 10^®(i?7S“^ (3) 


Stellar mass and age of starburst: The ratio 
g(H®)/Ly decreases rapidly with the mass of the most 
massive surviving O stars, i.e., the age of the starburst. 
Using the Starburst99 model of Leitherer et al. (1999) for 
metallicity 0.4 solar and normal Salpeter IMF, we deter¬ 
mine an age of 3.3 Myr. Once this is known, the Lyman 
continuum luminosity can be used to infer the total mass 
in the corresponding starburst; we obtain 400 Mq. The 
number of O stars involved is 1-2, so these numbers must 
be considered as averages. 

Mass and size of gas cloud: We can estimate the 
total mass of ionized hydrogen from (Osterbrock 1989) 

Mhii = Q(H°)mp/[neaB] = 420n/Qodi7 M©, (4) 

where nioo = ne/100cm“® and rrip denote the electron 
density and proton mass. Around early O stars, the H 
Stromgren radius is 


MHii(l-ky+) 


11/3 


(47r/3)nemp 


- pc, (5) 


mil = 
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where y~^ ~ 0.1 is the fraction by number of ionized he¬ 
lium. This is indeed much smaller than the photometric 
spatial resolution (about 0".7 = 57di7pc). The column 
density is 

IVhii = Mmi/ (ttchii) = 1.3 x cm“^. (6) 

From the inferred extinction and metallicity we may es¬ 
timate the intervening hydrogen column density from the 
local interstellar medium relation (Bohlin, Savage & Drake 
1978), using 

1V(H) ~ 5.9x10^^ cm“^mag“^EB_yZ0/Z ~ 6x10^^ cm“^. 

If the neutral hydrogen is in a much denser shell than the 
HII, then we can use the HII radius also to estimate its 
total mass 

Mh = Trr^ii X 2iV(H)mp = 3600 M©. (8) 

The star formation efficiency would then be around 10%. 

4. DISCUSSION 

The compact Hll-region is located about 3.'4 (17 kpc 
projected distance) north and 0.'9 (4.4 kpc) west of NGC 
4388, almost perpendicular to the disk plane of this galaxy 
and 45° away from the nearest part of the very extended 
emission-line region (VEELR) discussed by Yoshida et al. 
(2002). The radial velocity inferred from the emission lines 
is 2670kms“^, whereas the galaxy has Vr = 2520kms“^. 
The near-coincidence of these numbers may indicate that 
the HII region is, or once was, part of the NGC 4388 sys¬ 
tem. From the large radial velocity relative to the Virgo 
center {vr ~ l.Scryirgo) and from the Tully-Fisher distance 
of NGC 4388 (Yasuda et al. 1997), both are probably 
falling through the cluster core. 

This HII region is powered by a small stellar associa¬ 
tion or star cluster, with an estimated mass of ~ 400 Mq 
and age of ~ 3Myr for a normal IMF. If it is a young 
cluster, it must have a radius smaller than that of the HII 
region (~ 3.5pc). Clusters with these parameters have 
short relaxation times and dissolve by internal dynamical 
processes; within a few 10® yr the stars would be added to 
the diffuse stellar population nearby. 

The formation of such objects is thus possible far out¬ 
side the main star forming regions in galaxies. Perhaps the 
most plausible explanation for the observed position and 
velocity of the HII region relative to NGC 4388 is that it 
is already unbound and moving on a different orbit in the 
cluster potential. In this case we would be seeing intraclus¬ 
ter star formation. If on the other hand its true distance 
from NGC 4388 is comparable to the projected distance 
so that it is still bound, we would have discovered a small 
star-forming knot in the far outer halo of NGC 4388 far 
from any other star formation activity. This would have a 
velocity in the frame of the galaxy comparable to the cir¬ 
cular velocity (~ 200kms“^) of NGC 4388. In fact, in the 
Subaru field that contained our ICHH region, Arnaboldi 
et al. (2002) found a small sample of similar candidate 
objects, among which 7 candidates are located in the out¬ 
skirts of M86, at a distance of 10-15 kpc in a disturbed 
region which probably also contains diffuse Ha emission, 
and 1 such object in the outer parts of M84. 

What has triggered the recent onset of star formation in 
this remote cloud? Notice that with its radial velocity, it 


will have moved ~ 4 kpc relative to the Virgo cluster, and 
~ 500 pc relative to NGC 4388, in the lifetime of the mas¬ 
sive stars. Thus these stars will essentially have formed in 
situ. A possible model could be that the cloud was com¬ 
pressed after entering the hot intracluster medium (ICM). 
While the typical pressure of the ICM {ueT ~ lO'^ K/ cm®) 
is far smaller than the internal pressure inferred from typ¬ 
ical HII region parameters(ne ~ 100/ cm®; T ~ 10"^ K), it 
could be comparable to that of the surrounding neutral or 
molecular cloud traced by the absorption. The ionized re¬ 
gion may be reexpanding into the surrounding cloud and 
into the ICM - the mass of the cluster is not large enough 
to bind the HII region. 

Reanalysing the Subaru field photometry, we have also 
detected a number of extended emitters in both [OHI] and 
Ha, with colors similar to the confirmed compact HII re¬ 
gion. These are peaks in more extended emissions, and lie 
at distances 11-33 kpc along the direction of the VEELR 
identified by Yoshida et al. (2002) towards the NE of 
NGC 4388. These objects have relatively low excitation 
([OHI]/Ha 0.6, in agreement with Fig. 6 of Yoshida et 
al. 2002), and show continuum emission near the limiting 
magnitude of the combined (V-l-R) image (Arnaboldi et al. 
2002). Their half-light radii are similar in the combined 
continuum and in the [OHI] and Ha images, which fur¬ 
ther supports the notion that these distant sources have 
underlying continuum emission, and are thus ionized by 
OB stars rather than by the nucleus of NGC 4388. From 
the discussion of Yoshida et al. (2002) it is likely that this 
star formation is induced by the jet or ionization cone. 

Yoshida et al. (2002) proposed that the VEELR ob¬ 
served to the NE of NGC 4388 could be the debris of a 
past interaction. Although our ICHH region is 2.6 arcmin 
away from the nearest part of the VEELR, it could also be 
associated with such tidal debris. The gas in this ICHH 
region seems unlikely to come from ram pressure stripping 
of the gas in NGC 4388, because the peculiar velocity of 
the ICHH region relative to the systemic velocity of the 
Virgo cluster is even larger than for NGC 4388 itself. (We 
thank B. Moore for this comment.) 

The total intracluster star formation rate (ISFR) es¬ 
timated from the candidates in this field is small, how¬ 
ever. Counting the 8 ICHH candidates near M86 and 
M84 plus 3 intracluster candidates, including that stud¬ 
ied here, and the 6 best extended candidates, gives a to¬ 
tal of 17 star-forming HII regions in the Subaru field. If 
the density of these objects is typical for the Virgo clus¬ 
ter core, these would correspond to ~ 10® such objects 
throughout Virgo. If we adopt similar parameters to those 
inferred for the HII region here (conservatively), we ob¬ 
tain from these candidates a total ISFR of ~ 10 times 
~ 400 Mq per 3Myr, in a surveyed area of 918 arcmin^, 
i.e., ~ 10“® M0arcmin“^ yr“®. For comparison, the intra¬ 
cluster luminosity inferred from planetary nebulae in the 
Subaru field is ~ lO^LB^0arcmin“^. However, it appears 
not impossible that at higher redshifts, when the environ¬ 
ment of infalling galaxies was more gas-rich, formation of 
stars directly from this intracluster gas could have been an 
important part of the origin of intracluster stars in Virgo. 

The existence of compact HII regions in Virgo, beyond 
its intrinsic interest, is relevant for a few other issues as 
well. Firstly, the massive stars ionizing the gas explode as 
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type II supernovae and enrich the Virgo ICM with met¬ 
als. This process adds to the main metal content of the 
ICM, which is believed to have occured at high redshift 
when most of the stars in the cluster ellipticals and bulges 
formed (Renzini 1999). From the ISFR inferred above, the 
present expected supernova rate is ~ 10“®arcmin~ yr“^. 
SN II from isolated star formation could provide the en¬ 
richment inferred for the Lya clouds in Virgo (Tripp et 
al. 2002). Because our ICHII region contains only 1-2 O 
stars, observations of the enrichment from such objects 
in isolated regions might give constraints on the element 
yields of individual supernovae. Newly formed stars from 
this material could then have a range of abundance ratios 
similar to old Galactic halo stars (Argast et al. 2000). 

Secondly, compact HII regions may affect distance deter¬ 
minations with the planetary nebula luminosity function 
(PNLF, Jacoby et al. 1992). The [OIII] luminosity of this 
ICHII region places it at the bright end of the luminosity 
function inferred from [OIII] emission line surveys. In a 
sufficiently deep offband control image its continuum light 
would be visible, leading to removal from the PN sample. 
However, this does require a deep image, and moreover 


there are a few HII candidates in the Subaru field for which 
the continuum is not detected. A few unrecognized such 
objects per galaxy or intracluster field could cause the dis¬ 
tance to these fields to be underestimated. It is unclear 
whether this may account for the lower average distance 
inferred by the PNLF compared to surface brightness fluc¬ 
tuations ( Ciardullo et al. 2002). 

Thirdly, if compact ICHII regions exist in galaxies gen¬ 
erally, they could be the birth places of distant B stars 
in the Galactic halo, some of which are too far from the 
disk to have been ejected from there into the halo (Gon- 
lon et al. 1992). Most likely, their birth places would be 
distant high velocity clouds. At this time, however, there 
is no evidence for star formation in galactic high-velocity 
clouds. 


We thank R. Scarpa for efficient help at UT4, J. Alcala 
for independently checking line fluxes, and A. Capetti for 
a helpful discussion. OG thanks the Swiss Nationalfonds 
for support under grant 20-64856.01. This research has 
made use of the NASA extragalactic database. 
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